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Abstract. Explicit factorized formulas for the matrix elements (form- factors) of the 
spin operators a x and er y between the eigenvectors of the Hamiltonian of the finite 
quantum periodic XY-chain in a transverse field were derived. The derivation is based 
on the relations between three models: the model of quantum XY-chain, Ising model 
on 2D lattice and N = 2 Baxter-Bazhanov-Stroganov r^-model. Due to these 
relations we transfer the formulas for the form-factors of the latter model recently 
obtained by the use of separation of variables method to the model of quantum XY- 
chain. Hopefully, the formulas for the form-factors will help in analysis of multipoint 
dynamic correlation functions at a finite temperature. As an example, we re-derive 
the asymptotics of two-point correlation function in the disordered phase without the 
use of the Toeplitz determinants and the Wicncr-Hopf factorization method. 
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1. Introduction 

The quantum XY-chain is one of the simplest models which is rich enough from the point 
of view of physics and at the same time permits strict mathematical analysis. The study 
of this model was started in pQ where it was rewritten in terms of fermionic operators 
by means of the Jordan- Wigner transformation. Now this relation is a standard mean 
to study different properties (the spectrum of the Hamiltonian [H [2], the correlation 
functions 0, BJ El El [7] , the emptiness formation probability [B], the entanglement entropy 
0, [10], [HI [121 EH]), quantum quenches [H] in XY-chain. Although the Hamiltonian of 
the model is equivalent to the Hamiltonian of a free fermionic system, the spin operators 
cr x and a 7 are expressed in terms of fermionic operators in a non-local way. Thus the 
study of correlation functions of such operators is a non-trivial problem. For example, 
the correlation function (er x er x ) can be written through Toeplitz determinant of size d 
and the derivation of the asymptotics d — > oo requires [H] the use of Szego theorem and 
the Wiener-Hopf factorization method. 

In this paper we propose an alternative way to study correlation functions of the 
XY-model: we derive the formulas for the matrix elements of spin operators <r x and 
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cr y between the eigenvectors of the Hamiltonian of the finite quantum XY-chain in a 
transverse field. These formulas allow to obtain at least formal expression for multipoint 
dynamic correlation functions at a finite temperature. For this aim, it is enough to 
insert the resolution of identity operator as a sum of projectors to the eigenspaces of 
Hamiltonian. Hopefully the correlation functions in terms of these sums will be more 
easily analyzed. As an application of the formulas for form-factors, we re-derive the 
asymptotics of correlation function (c^a^) at d — > oo. 

The idea of derivation of form-factors of the quantum finite XY-chain is to use the 
relations between three models: the model of quantum XY-chain, the Ising model on 
2D lattice and N = 2 Baxter-Bazhanov-Stroganov (BBS) model [T5l [16] . The relation 
between the first and the second model was observed in [17] (the relation ffl6|) in this 
paper between the energies of fermionic excitations of these two models seems to be 
new), the relation between the second and the third model was found in (T8J [19]. The 
parameters of the models are (h,x), (K x ,K y ) and (a, b), respectively. Due to these 
relations we transfer the formulas for the form-factors of N = 2 BBS model, recently 
obtained [201 121] by the use of separation of variables method, to the model of quantum 
XY-chain. The main formulas are flHTjl . (152) together with fl2TT]), (1291), (jlBjl . 

The separation of variables method for the quantum integrable systems (with basic 
example being the Toda chain) was introduced by Sklyanin [22] and further developed 
by Kharchev and Lebedev [23]. In [19], this method was adapted for BBS model (a Z_/v- 
symmetric quantum spin system) to obtain the eigenvectors of this model. At N = 2 and 
special values of parameters the BBS model reduces [HJ EHJ El] to the Ising model. The 
eigenvectors of transfer-matrix of Ising model obtained by the separation of variables 
method allowed [20lEl] to prove the conjectural formula j25j|26] for the matrix elements 
of spin operator for finite Ising model. This derivation had provided first proof of the 
formula. A summarizing overview of the results on separation of variables for BBS 
model is given in [21] . It is interesting that factorized formulas for the matrix elements 
of spin operators exist also for superintegrable Z^-symmetric chiral Potts quantum chain 
[271128]. 

In Sect. 2 we remind the definition of the finite quantum XY-chain in a transverse 
field, its phase diagram, eigenvalues of the Hamiltonian and give general comments on 
the matrix elements of spin operators between the eigenvectors of the Hamiltonian. 
Sect. 3 is devoted to the description of relations between three models: the model 
of quantum XY-chain, the Ising model on 2D lattice and N = 2 Baxter-Bazhanov- 
Stroganov model. Using these relations, in Sect. 4 we derive formulas for the matrix 
elements (form-factors) of the spin operators o~ x and cr y between the eigenvectors of the 
Hamiltonian of the finite quantum XY-chain. In Sect. 5, these formulas are rewritten 
to the case of chain of infinite length. In Sect. 6, as an application of the formulas for 
form-factors, we re-derive the asymptotics [3] of correlation function (c^a^) at d — > oo 
without the use of the Toeplitz determinants and Wiener-Hopf factorization method. 
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2. Definition of the finite quantum XY-chain in a transverse field 

2.1. The Hamiltonian and phase diagram 

The Hamiltonian of the XY-chain of length n in a transverse field h is [H [2] 

n = ~\ £ (^^x+i + ^^x+i + * *s) , (i) 
k=i ^ ^ 

where a\ are Pauli matrices, x is the anisotropy. In the case x = we get XX-chain 
(isotropic case). The value x = 1 corresponds to the quantum Ising chain in a transverse 
field. In what follows we restrict ourselves to the case x > 0, h > 0. Other signs of x 
and h can be obtained using automorphisms of the algebra of Pauli matrices. Also we 
will suppose the periodic boundary condition a l k = o~ k+n . In [29] it is shown that the 
formulas for the matrix elements of spin operators obtained in Sect. H] are also applicable 
for the antiperiodic boundary condition o~ k+n = —o~ k , o~ y k+n = —a\ and o~ k+n = o\. 

Now about the values of h. Due to the relation of XY chain with 2D Ising model, 
which will be discussed in the next section, the coupling constant h plays the role 
of a temperature-like variable. The value h > 1 corresponds to the paramagnetic 
(disordered) phase. The value < h < 1 corresponds to the ferromagnetic (ordered) 
phase. If < h < (1 — x 2 ) 1 / 2 , it is oscillatory region (because of oscillatory behavior 
of two-point correlation function). Another peculiarity related to this region is the 
following. At fixed x, < x < 1, in the region where h > (1 - x 2 ) 1 / 2 the NS- 
vacuum energy is lower than R- vacuum energy (asymptotically, at n — > oo, they become 
coinciding). In the region < h < (1 - x 2 ) 1 / 2 there are intersections at special values 
of h of these vacuum levels even at finite n. The number of these intersections grows 
with n. For a detailed analysis of the oscillatory region see [31 [30] . 

In this paper we derive the formulas for the matrix elements in paramagnetic phase 
(h > 1, < x < 1) and add comments on the modification of the formulas for other 
values of parameters. 

2.2. Eigenvalues and eigenvectors of the Hamiltonian of XY-chain 

Using Jordan-Wigner and Bogoliubov transformations the Hamiltonian "H of the XY- 
chain can be rewritten as the Hamiltonian of the system of free fermions and diagonalized 
[U|2]. The relation between energies e(q) and momenta q of the fermionic excitations is 

e(q) = ((h - cosg) 2 + x 2 sin 2 g) 1/2 , g^0,7r, (2) 

e(0) = h-l, e{n) = h + l. 

The Hamiltonian % commutes with the operator V = o\o\ • • -a^. Since V 2 = 1, 
the eigenvectors are separated to two sectors with respect to the eigenvalue of V. 
Below the sign +/— in front of e(q) in the expression for energies S corresponds to 
the absence/presence of the fermionic excitation with the momentum q. Each such 
excitation carries the energy e(q). 
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• NS-sector: V — > +1, the fermionic excitations have "half-integer" quasimomenta 

q € NS = [— (j + 1/2) , j e Z n \ =► 5 = -1 Yl ^ ■ 

This sector includes the states only with an even number of excitations. 

• R-sector: V — > — 1, the fermionic excitations have "integer" quasimomenta 

?£R=ft', iezJ E = -~J2±<q)- (3) 

In the paramagnetic phase this sector includes the states only with an odd number 
of excitations. In the ferromagentic phase (0 < h < 1) it is natural to re-define 
the energy of zero-momentum excitation as e(0) = 1 — h to be positive. From the 
formula ([3]) for energy S, this change of the sign of e(0) in the ferromagnetic phase 
leads to a formal change between absence/presence of zero- momentum excitation in 
the labelling of eigenstates. Thus although the analytical expressions for energies 
8 in terms of h and x are the same in both phases, because of the redefinition of 
5(0) in the case of < h < 1 the number of excitations in the ferromagnetic phase 
is even. 

We will denote the eigenstates |$)modci by the set of values of the excited quasi- 
momenta $ = {qi, q 2 , ■ ■ ■ , qi}, the label of the sector a = NS or R and the label of the 
model. For example, the state in R-sector of the quantum XY-chain with n = 3 sites in 
the paramagnetic phase with all possible quasi-momenta $ = {0, 27r/3, — 27r/3} excited 
is |0, 27r/3, — 27r/3)^ y . It has the energy 

£ = --(l-h- £(2tt/3) - e(-27r/3)) . 

The same formula for the energy in the ferromagnetic phase corresponds to the state 
|27r/3,-27r/3)* y . 

2. 3. Matrix elements of spin operators 

Formally in order to calculate any correlation function for XY-chain it is sufficient to 
find the matrix elements of spin operators a£, a\ and o~\ between the eigenstates of 
Hamiltonian %. 

• Matrix elements of a\. 

The operator o~\ commutes with V = o\a\ ■ ■ ■ a^. Therefore the action of a\ does 
not change the sector. In fact the operator o\ can be presented as a bilinear 
combination of operators of creation and annihilation of the fermionic excitations. 
Hence the matrix elements of o~\ between eigenvectors of % can be calculated easily 
(most of them are 0). We will not consider such matrix elements in this paper. 

• Matrix elements of and o~\. 

The operators <r£ and o~\ anticommute with V = o~\o~\ ■ ■ ■ o~ z n . Therefore their action 
changes the sector. The operators 0% and a\ can not be presented in terms of 
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fermionic operators in a local way. All the matrix elements of them between the 
eigenvectors of H from different sectors are non-zero! 

The aim of this paper is to derive explicit factorized formula for the matrix elements 
of cr£ and a y k . The idea is to relate three models: the quantum XY-chain in a transverse 
field, the Ising model on 2D lattice and N = 2 BBS model. The relation between the 
first and the second model is based on the observation by M. Suzuki [17]. The relation 
between the second and the third is based on [18J. The latter relation together with the 
results on separation of variables for BBS model allowed [21] to prove the formulas for 
the matrix elements of spin operator of Ising model found by A. Bugrij and O. Lisovyy 
[251 126] . In this paper we transfer these results on the matrix elements to the case 
of XY-chain. The parameters of these three models are (h,x), (K x ,K y ) and (a, b), 
respectively. 



3. Relation between three models 



3.1. Relation between quantum XY-chain and the Ising model on a lattice 

The row-to-row transfer-matrix of the two-dimensional Ising model with parameters K x 
and K y can be chosen as 

/ n K* \ ( n \ ( n K* \ 

t XY := Tl' 2 T 2 Tl' 2 = exp PT -f a\ exp HT K x < exp HT -f a% , (4) 

\fc=i / \fc=i / \fc=i / 

where the spin configurations of the rows are chosen to be labeled by the eigenvectors 
of the operators <r£, the parameter K* is dual to K y , that is tanhi^ = exp(— 2K*), and 

Tl = ex P K l ^ J ' T2 = ex P k * a * °M ■ ( 5 ) 

In [I7J, M. Suzuki observed that if we choose K x and K* such that 



yj\ — x 2 1 

tanh 2K X = — r , cosh 2K* = - (6) 

h y x 

then the Hamiltonian ([1]) of XY-chain will commute with the transfer-matrix of the 2D 

Ising model 01]) and these two operators have a common set of eigenvectors. 

3.2. N = 2 BBS model and its relation to the Ising model 
To define N = 2 BBS model we use the following L-operatoJj] [T6l 131] 
r ( 1 +Xal X °Ua-bo-t) \ 

LkW = , (7) 

\c-Ua-bo-l) W+o%P J 

depending on parameters a, b and spectral parameter A. It satisfies the Yang-Baxter 
equation with (twisted) quantum trigonometric i?-matrix. In particular it means that 

| In comparison with [21] we interchanged a\ and a\ in the L-operator ([7|. It is just another 
representation of the Weyl algebra entering the definition of L-operator. 
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the eigenvectors of the transfer matrix t(A) = trLi(A) L 2 (X) ■ ■ ■ L n (X) built from such 
L-operators are independent of A. 

Fixing the spectral parameter to the value A = b/a, the L-operator ([7]) degenerates 



and the transfer matrix t(A) can be put into a nonsymmetric Ising form 

n n 

t(b/a) = + a% ■ b/a) ■ + o%o* +1 ■ ab) = (cosh cosh K* y )~ n , 

k=l k=l 

t(b/a) ~ T{T 2 = exp (ELi K a z k ) exp (ELi K* °l , (8) 
if we use periodic boundary condition o~ l n+k = a\ and identify 

e -2K y =tan hi^* = b/a, tanhK x = ab . (9) 

Thus at A = b/a we get the transfer-matrix of the Ising model. If we do not fix the 
spectral parameter to this special value, we shall talk of the "generalized Ising model" . 
However, transfer matrix eigenstates are independent of the choice of A. In [191 EDI EE] 
the eigenvectors for the nonsymmetric transfer matrix ([8]) and the matrix elements of 
o\ between these eigenvectors were derived using method of separation of varibales. 

Comparing fl6]) and we get the following simple relations for the parameters of 
XY- model and special BBS-model with L-operator (J7|): 

a 2 -b 2 , l + a 2 b 2 , , 

X= ^T^' h = lP + W (10) 

3.3. Relation between the energies of excitations for the Ising model and XY-chain 

In the previous subsections we have shown how the quantum XY-chain in a transverse 
field, the Ising model on 2D lattice and N = 2 BBS model are related. The parameters 
of the models are (h,x), (K x ,K y ) and (a,b), respectively. The relations between these 
pairs of parameters are given by (jSJ), © an d (ITU1) . 

Since the transfer- matrices txy = T^ 2 T 2 T^ 2 , ti s = T 2 ^ 2 TiT 2 ^ 2 , T\T 2 of 2D Ising 
model are related by similarity transformations, the enumeration of the eigenstates of all 
these transfer-matrices is the same as described in Sect. 2.2 for the quantum XY-chain. 
They will be denoted, respectively, by |$) X y> I^Sh \®) a i where $ = {pi,p 2 , ■ ■ ■ ,Pl} is 
the set of values of the excited quasi-momenta $ = {pi,p 2 , ■ ■ ■ ,Pl} an d a = NS or R is 
the label of the sector. Their eigenvalues e -7 * are the same: 

txY|$)xY = e- 7 *|$)xY, t*,|$)2 = e-»|$)£, T^]*)" = e~^ , (11) 
L 1 

(=1 pea 

cosMri ^+g^' -^os,,, (18) 
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t x = tanh K x , t y = tanh K y . 

The eigenvalues of the transfer matrix t(A) of the BBS model with L-operator ((7j) 
are proportional to Yl P (^ ^ s p) ( see formula (68) of [20]). 

( b A -2b 2 cosp+l \ 1/2 
Sp ~ U 4 -2a 2 cosp + ly' ' ( j 

where the sign +/— in the front of s p in the expression for the eigenvalues of t(A) 
corresponds to the absence/presence of the fermionic excitation with the momentum p. 
The momentum p runs over the same sets as in the case of the quantum XY-chain. Due 
to ([8]), f[TT]) and ([12"]) we have the relation between ^{p) and s p (see [21]): 

C 7G0 = ^L±l (15) 
as p — o 

and a relation between e{p) and 7(p): using (TlOl) . ( THj) and ([T5l) we get 



2ab(a 2 + b 2 ) , . VT^? ^ ( . R , 

■ e(p) = — ; s , , = e(p) • (16) 



(a 2 -6 2 )(l-a 2 6 2 ) v ^ x^TW 

The existence of the relation between 7Q9) and is surprising because the commu- 
tativity of the Hamiltonian (CQ) of the XY-chain and the transfer matrix ([4]) of the 2D 
Ising model does not imply a priori any relation between their eigenvalues. 



3.4- Uniformization of the dispersion relation [13]) 

We use a parametrization of the dispersion relation (TT5|) of the 2D Ising model in terms 
of elliptic function at/i>l,0<x<l which corresponds to the paramagnetic phase 
of Ising model. This parametrization is a modification of parametrization from [52] 
given for the ferromagnetic phase of Ising model and corresponding to < h < 1, 
(1 - h 2 fl 2 < x < 1 for XY-chain. 

We introduce the modulus of elliptic curve k -1 by 

k" 1 = sinh 2/^ sinh 2^ = x/Vx 2 + h 2 - 1 = (a 2 - b 2 )/(l - a 2 b 2 ) . 

In the paramagnetic phase we have < k _1 < 1. Complete elliptic integrals for k _1 and 
for the supplementary modulus are K = K(k~ 1 ) and K' = K((l — k -2 ) 1 / 2 ), respectively. 
We define real parameter a, < a < K'/2, by one of the equivalent relations 

l/sim^-fQ =ik/sn(2ia,k~ 1 ), l/wah2K v = -i sn(2ia, k" 1 ) . 

The following two elliptic functions 

X(u) = sn(u — ia, k )/sn(it + ia, k" 1 ) , z{u) = k" 1 sn(w — ia, k" 1 ) sn(w + ia, k _1 ) 
satisfy the relation 

sinh 2K X (z + z~ l ) + sinh 2K y (A + A -1 ) = 2 cosh 2K X cosh 2K y . (17) 
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To prove it we note that the left-hand side of (1171) is an elliptic function without poles 
(the poles at ±a and ±a + \K' are canceled) and therefore it is a constant. Thus it is 
sufficient to establish validity of ( TT7T) at u = 0. For this end we use A(0) = —1, 

,1 9, , u 1 - dn(2ia, k" 1 ) w , , 1 + dn(2ia, k" 1 ) 

z(0) = -k- 1 sn 2 ia, k" 1 = -k ; o . '- 1N ; , z' 1 = -k ) o . 1N ; 

w v ' 1 + cn(2ia, k" 1 ) K 1 1 - cn(2ia, k" 1 ) 

which follow from the formulas of Example 6, Sect. 22.21 of [33 J, and 

cosh2Ji x . = dn(2ia, k" 1 ) , cosh 2K y = i cn(2ia, k~ 1 )/sn(2ia, k~ x ) . 

The relation ( TTTj) coincides with the dispersion relation (Tl3|) if one identifies 
z(u) = e~ ip and X(u) = e~"'( p \ The parameter u on the elliptic curve is an analogue 
of rapidity. Now if p runs from — 7r to 7r then u runs along the segment from iK' /2 to 
2K + iK'/2. 

There is another dispersion relation corresponding to the evolution in the transverse 
direction on the Ising lattice: 

(U + 1- 1 )^ + 1- 1 ) t~ l - 1 _ 
cosh7(p) = — j- r-= — cosp. (18) 

It is uniformized by \(u) = e ip and z(u) = e~^^\ Now if p runs from — n to n then u 
runs along the segment from to 2K. 
From (fTUj) . we have 

2 h - Vh 2 + x 2 - 1 „ /i - Vfr 2 + x 2 - 1 

a = , = — . (19) 

1 — x 1 + x 



The points with z = a ±2 and z = b ±2 are the branching points of the spectral curve ((Tj 
considered as \{z). The parameters a 2 , b 2 correspond respectively to A2, Aj" 1 of [3] and 
to «r\ of [32]. We have also a 2 = e~^°\ b 2 = e~^ n \ 



4. Formulas for the matrix elements of spin operators 

In this section we will derive the formulas for the matrix elements of spin operators for 
quantum XY-chain of finite length. The derivation for the basic region of parameters 
/i>l,0<x<lis given in Sect. 14.11 The formulas for other values of parameters can 
be obtained by analytic continuation. The details of the continuation are given in the 
following subsections. 



4-1. Paramagnetic phase: h>l, 0<x<l 

We use the Bugrij-Lisovyy formula ((40) of [2B]) for the matrix element of spin operator 
between the eigenvectors |$ )is = \Qi, <?2, • • • , qi<)f s S and |$i)i s = \pup 2 , ■ ■ ■ ,Pl)u of tne 
symmetric transfer matrix ti s = T^ 2 TiT^ 2 for the finite 2D Ising model (the states are 
labeled by the momenta of excited fermions as it is explained in Sects. 12.21 and 3.3): 

^$0,$! = I (qi, q-zi ■ ■ ■ , Qk I o"m I Pi,P2> ■■■iPl )f s \ 2 
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K T7 NS ~\-nh Tfa)+7(g) L ttR • i l{pi)+l(p) /, - -1 \ {K-Lf /2 

A W<?fc 2 tt 1 Ip^pi blim 2 ( l y~ l y 



^ T II „nR„ ; „v,7(ft)+7W n 



„ Hp sinh 2toM 11 „ J]NS ginh 7fa)+7(<?) ^ _ t-i 



X 



TT sm ~ 2~ TT sm ~2~ TT smil 2 f20 N 
11 ^^TW+lM II^tW+iW 11 sin 2 Mi ' 1 ; 

k<k> S1Iin 2 l<l> S1Iin 2 ^<k<K 0111 2 



£ 4 = k 2 - 1 = (sinh 2K X sinh 2i^)" 2 - 1 = — = V , = ; \ 2 , 2 ' , (21) 



a 2 - 6 2 ) 5 



n NS ginh tWMI ^ sinh 7M±7M 



ta; = tanh = ab , t y = tanh Kj, 



a + 6 



t„ - V 1 4a 2 6 2 1 - x 2 



ic-t- 1 (a2-62 )(1 _ a 2 6 2) Xv / X 2 + h 2 _ 1 ' (23) 

where we used ([6]), (J9J) and (fTUl) to write equivalent expressions in terms of different 
parameters. 

Note that the Bugrij-Lisovyy formula ( [20]) is given for the normalized eigenvectors 
|$)i s of the transfer matrix tj s = T^ 2 TiT^ 2 which differs from the transfer matrix (j4j). 
The eigenvectors |$)xy of the Hamiltonian ([I]) of XY-chain and (J4j) are in one-to-one 
correspondence with the eigenvectors |$)i s and the eigenvectors |<&) of BBS model. All 
the eigenvectors with the same $ (the same set of excited fermion excitations) are related 
by similarity transformations. 

In [191 120] the left and right eigenvectors (<&| and |$) of T{T 2 were found. They 

1 /2 

are related to |$)i s by the action of operator T 2 and its inverse. Since these operators 
commute with <r^, we have the relation 





" m 








($0 


$0>($l| 


$1> 



(24) 



expressing the matrix elements of 2D Ising model in terms of the matrix elements of 
BBS model found in [21] and used to prove (12"0|) . 

Here we want to use the matrix elements of BBS model to derive the matrix elements 
of spin operators between the eigenstates of XY quantum chain. Since the Hamiltonian 
of XY-chain and the transfer matrix ti s are Hermitian matrices, there is a natural way to 
relate the left and right eigenvectors and to normalize them. In the case of nonsymmetric 
transfer matrix T1T2 arising in BBS model, the left and right eigenvectors ($| and 
|$) are unrelated but we will relate them to normalized eigenvectors of XY-chain by 
a\ ■ xy($ I = ($|Ti /2 and a% ■ |$) XY = T[ 1/2 \<S>). We fix a% = ||Tr 1/2 |$}|| > 0, then 
the coefficient a\ is determined from the requirement xy( < ^ ) | = |^ ) )xy- We have also 
ot\a% = ($|$). Since ($|$) is real at real a and b, both a\ and a$ are real too. The 
matrix elements of BBS model are related to the matrix elements of XY-chain by 

(*oKI*i> = « • xY($o|Tr 1/2 ^T 1 1/2 |$ 1 } XY (25) 
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where the last relation follows from the facts that |$)xy is eigenvector of T^ 2 T 2 T^ 2 
with the eigenvalue e~ 7 * (see (1TTT) ) and T 2 commutes with cr^. Complex conjugation of 
f )26|) together with (125]) with interchanged $o an d ^1 gi ye 

L R 



<*oKI*i> = • e^o-7* lxY ($ | T vv^ i Tr 1/ 1$ 1 )xY • (27) 



From the other side, the formulas (56) and (57) of [2T] give the following factorized 
presentations for the matrix element between eigenstates of the transfer matrix TiT 2 
(see footnote on p. [5]): 

($o|<J$i) = h(b)f 2 (b 2 ) , ($iK|$o> = h(-b)'W), 

where fi(b) is a real (we suppose that a and b are real) and f2{b 2 ) is a complex but 
invariant with respect to b — > —b. Therefore 

($oKI$i) Mb) _ r r9R , 
It is easy to prove using explicit formulas for fi(b) from [21] that 

for |$ ) = \Qi, 92) • • • j <?_ft:) NS and |$i) = |pi,p 2 > • • • ,Pl) R - Let us consider, for example, 
the case of odd n (the length of the chain) and o~ = (a q = 0/a q = 1 corresponds 
to absence/presence of fermion excitation with momentum q). From Eq. (56) and the 
discussion in Section 6.1 of [21] we can choose 

h{b) = ((-ir -ab) l[((-l) k b + as nk/n ) ]J((-l) k b - as« k/n ) , 

ket> kev 

where s p is given by ( 1T%|) and the set V (resp. V) consists of such k from {1, 2, . . . , n— 1} 
for which the fermions with both momenta ±irk/n are excited (resp. not excited) in the 
states |$ ) an d Using (fTB]) we get 



C, 



a(6) _(-ir -06 ru, P ^ 7(p)/2 nf=i, 9fe ^- 7(9fe) 



e ' 



Thus it remains to verify that the first fraction also fits ( 1291) as contribution of the 
momenta and tt. For this end we need just to take into account (fl4|) . ([TBI) and 

a 2 -l' w a 2 + l' "l + afe' 

The correct sign of the latter formula can be fixed from the limit 6 = (the quantum 
Ising chain limit) and then taking limit a — > (it corresponds to the limit of strong 
external field h — > 00). All the other three cases of odd (even) n and cr = °~n (°o ^w) 
can be analyzed similarly. It proves (1291) . 
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Taking into account ([M]) and (1251) we get 



(($ |$o)($i|$i)) 1/2 

where ^$0,$! is a phase related to a particular normalization of eigenvectors. To relate 
these matrix elements to the matrix elements of XY-chain we observe that (127]) and 
(J2SD imply a^o^ = a^a^. Since a^a^a^a^ = ($ |*o)(*i|*i) we obtain 
c^R = (($ |$ )($ 1 |$ 1 )) 1 / 2 . Thus from ([25]), fl2E), (ESJ) and (ED]) we derive 

x Y ($o|Tr 1/2 <T 1 1/2 |$ 1 ) XY 

= xy($o| «coshiT* - i< sinhK;) I^xy = e i5 *o,*i (C $0)$1 H$ 0)$1 ) 1/2 , 
x Y ($o|T 1 1/2 <Tr 1/2 |$ 1 ) XY 

= xy($o| «cosh^ + ia y m smhK;) I^xy = e w *°^ (C^ S 0O , $1 ) 1/2 . 

Finally taking appropriate linear combinations of these two formulas we get the main 
result of the paper: the matrix elements of spin operators between the eigenvectors 

|$o)xy = \qi, ?2, • • • , qi<)xY fr° m tne NS-sector and |$i)xy = \pi,P2, ■ ■ ■ ,£>l)xy from tlie 
R-sector of the Hamiltonian ([I]) of XY-chain are 

|xY($okml $ l)xY| 2 = 2( - 1 ~ ^ (cJ^Si + O^,*! ) S *o,*i 

cosh ^* ,*i j (31) 



1 + x 2 



|xY($o|^ml*l)xY| 2 - 2^ _ x (C*(, 2 *i ~ ^o,*i ) S *o,*i 



smh =#0,$! , (32) 



1 - x 2 

where S$ 0i $ i; C$ 0> $ 1 and 7$ are given by (120]) . (129]) . fTT2]) and (TT6]) . We also used 



^.2. Ferromagnetic phase: < x < 1, a/1 — x 2 < /i < 1 

Let us give a general comment on the continuation of the formulas from the region h > 1 
to region < h < 1. Formally, all the formulas for matrix elements of spin operators are 
correct for the paramagnetic phase where h > 1 and for the ferromagnetic phase where 
< h < 1. But for the case < /i < 1 it is natural to change the sign of e(0) (and also of 
7(0)) of zero-momentum excitation to be positive: e(0) = 1 — h. From ([3]), this change 
of the sign of e(0) (and of 7(0)) in the ferromagnetic phase leads to a formal change 
between absence/presence of zero-momentum excitation in the labelling of eigenstates. 
Therefore the number of the excitations in each sector (NS and R) becomes even. Direct 
calculation shows that the change of the sign of 7(0) in (13~TT) and ( 1321) can be absorbed 
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to obtain formally the same formulas (1t?Tj) and (I32p but with new 7(0), even L (the 
number of the excitations in R-sector) and new £ = (1 — k 2 ) 1 / 4 = ((1 — h 2 )/ 'x 2 ) 1 / 4 . 

The explicit formulas for the matrix elements for the region of parameters < x < 
1, y/1 — x 2 < /i < 1 are given 



4-3. Region x > 1 

From the relation ( I16p between the energies of XY-chain and Ising model excitations it 
follows that the energies i(jp) of the Ising model excitations become complex and it is 
useful to introduce ^(p) = vy(p) such that 



sin 7(2?) = l^ =£(p) . 

K ' xVx 2 + h 2 - 1 ^ ' 

Here should be chosen to be monotonically increasing function of p when p runs 

from to 7r. If h > x 2 — 1, the energy e(p) is monotonically increasing function of 

p with minimum at p = and maximum at p = tt. In this case < 7Q0) < tt/2. If 

< h < x 2 — 1, the energy e(p) is non-monotonical function having additional extremum 

(maximum) at p = p c , cosp c = h/(l — x 2 ), 7(p c ) = tt/2. In this case < 7Q0) < 7r/2 

for < p < p c and ^(p) > tt/2 for p c < p < tt. 

We continue analytically the formulas (l3~Tj) and (|32|) and write them in terms of 

7(p). All the changes in the final formulas are the following: 

1 , A 1 X i*\ 12 2x 2 7$ — 7$ ^ 
\xy(%K\<5>i)xy\ 2 = — — r cos 2 7 , 

X + 1 z 

|xy($oKI$i)xy| = rsm ^$o,1>i> 

x — 1 z 

where S$ <1 , 1 is given by (1211 with substitutions 

. , 7(P) + 7(?) , . T(p) + 7(9) ^ ~ ^ , x2 - 1 
sinn )• sin , — — > 



2 2 ' t x -t" 1 xVx 2 + h 2 - 1 

and 7$ is given by f fT2|) with ^{p) — > j(p). Also, for < h < 1, one has to take into 
account the modifications related to the zero mode described in Sect. 14.21 



4-4- Oscillatory region 0<x<l, < h < \/\ — x 2 

Similarly to the region x > 1, the energies 7(p) of the Ising model excitations become 
complex and it is useful to rewrite the matrix elements of spin operators in terms of 
lip) = lip) +i7r/2: 



cosh7(») = ejp) . 

Here ^(p) should be chosen to be monotonically increasing function of p when p runs 
from to tt. If 1 — x 2 < h < \/l — x 2 , the energy e{p) is monotonically increasing 
function of p with minimum at p = and maximum at p = tt. In this case 7Q9) > 0. If 
< h < 1 — x 2 , the energy e(p) is non-monotonical function having additional extremum 
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(minimum) at p = p c , cosp c = hj (1 — x 2 ), j(p c ) = 0. In this case j(p) < for < p < p c 
and 7(p) > for p c < p < n \ We continue analytically the formulas ( 13T|) and fl32|) and 
write them in terms of 7<jo). All the changes in the final formulas are the following: 

|xy<$o|^|^i)xy| = TTTq— — V e 2 +(-!) 2 e 



2(1 + x 



2 



Ixy^oKJ^xyI - tj-s \\ e 2 2 e 



2(1 



where H$ 0)$1 is given by ([20]) with £ = (1 - k 2 ) 1 / 4 = ((1 - /i^/x 2 ) 1 / 4 (since it is the 
ferromagnetic phase) and with the repacements 

sinh 7(g) + 7(g) ^ cosh 7(p) + 7(g) Vl£^ i-- 2 



t* - £ 1 Wi - x 2 - h 2 ' 

7$ is given by ( TT2|) with 7Q9) — > 7(p). Also one has to take into account the modifications 
related to the zero mode described in Sect. 14.21 



4-5. Other values of parameters 

4-5.1. Quantum Ising chain: x = 1. In the case of the quantum Ising chain (x = 1) 
the formula for the matrix elements of spin operator <r^ can be derived by a limiting 
procedure. The final formula was derived in [211 El] and it is expressed in terms of the 
energies of excitations e(q). In the case of general XY-chain we were not able to find an 
analogous formula for the matrix elements in terms of e(q). 

4-5.2. Boundary of oscillator region: x 2 + /i 2 = 1. One of the peculiarities of XY-chain 
when the parameters belong to the curve x 2 + h 2 = 1 is that the ground states |)x Y an d 
|)xy are degenerate and each of them can be presented as a sum of two pure tensors 
|35j . In [3], for such values of parameters, the two-point correlation function was found 
explicitly. Here we give some comments on the matrix elements of spin operators. They 
can be derived from the general formulas for the region < x < 1, yl — x 2 < h < 1 by 
a limiting procedure. From ( fl~6l) . denoting £ = a/x 2 + h 2 — 1 we get in the limit £ — > 

2h 

e7(p) = : n £ (p) ' £ 0°) = 1 ~~ h cos P ' 
Cv 1 - h 2 

- , lip) + 7(g) e^)+^))/2 

sinh = , £ = 1 , 6r = 1 ■ 

Using these formulas it is easy to take the limit ( — >■ in the general formulas for the 
matrix elements. We get, in particular, that the matrix elements are non-zero if and 
only if K = L or K - L = ±2. 
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5. Asymptotics of form-factors in the limit of infinite chain 

In this section we analyze the asymptotics of different parts of form-factors in the limit 
of infinite length (n — > oo) of XY-chain. For this end it is helpful to use the following 
integral representations for different parts of form-factors at finite n [26]. For 

and £t (see ( l22|) ) we have 

1 r j„ W 7(P) 



1 Z" 71 

A -1 = — dp logcoth 



1 T cosp-e n7(p) 
77(g) — — dp — — logcoth 



7r J cosh 7(g) — cosp 2 

n 2 C /" r dp dq*f'(pW(q) , sin((p + g)/2) 

log 



2vr 2 Jo </o sinh 717(79) sinh 777(7?) 



sin((p- g)/2) 



Let us show that A 1 — ^ if n — > 00. In fact we have the following asymptotics 
A" 1 -*- f dpe~ n ^ = - / 7r d pe -^(7(o)+7"(o) P V2+-) 



e 



T ./o 

717(0) 



7r y n7T7"(0) 

where we used the fact that 7(p) > (there is a gap in the spectrum for the non-critical 
parameters). Similarly we get 77(g) — > 0, £r — > 1 at n — > 00. For the derivation of these 
formulas together with the more precise asymptotics at n — > 00, see [36]. Another way 
to get the asymptotics in the limit of infinite length of the chain is given in [37] . 

In the limit of infinite XY-chain the formulas ( 13TT) . ( 1321) for the matrix elements of 
spin operators between the eigenstates |$o)xy = \Qi,Q2, ■ ■ ■ ,Qk)xy from the NS-sector 
and |$i)xy = \Pi,P2, ■ ■ ■ iPl)xy from the R-sector, and (1201) become 

I (ft \ n x i a \ ,2 _ ~ 2X pn ,u2 Ef=l7(gfc)-Ef = l7fa) ^ 
|XY\*0F m |*l/XY| ""Win COSjl o ' y 64 ) 

1 + x Z 

I |2_^ 2x •_ h 2 Ef=l 7(gfe) ~ Ef=l 7fa) r ,cN 

|XY\*0Fml*l/XY| — ^* ,*i 1 smn o j l d£) J 

1 — X z 



1 A 1 ^-c 1 ^ 72 



^ n n Rinti ^/(<77 n 



L = A nsinh7(g fc ) Al nsinh7(p { ) \* x — t" 1 



II • 1 2 7(q fc )+7(^) II .i2 7(P»)+7(P|/) II sin 2 ' ( 36 ) 

fc<fc' SmJ1 2 KV SmJ1 2 i<*<^ 2 



where the relation ( ITB"]) between 7(g) and the energy of excitation e(q) in XY-chain and 
relation (12"31 are also to be used. 



Form-factors of the finite quantum XY-chain 



15 



In the ferromagnetic phase (0 < h < 1) after an appropriate modification as 
explained in Sect. I4.2[ the thermodynamic limit formulas f )M|) and fl35|) at = L = 
allow to re-obtain [37] formulas for the spontaneous magnetization found in [T7] . Indeed, 
in the thermodynamic limit n — > oo, the energies of vacuum states |$o)xy = |)xy an d 
|$i)xy = |)xy °f different sectors asymptotically coincide giving the degeneration of the 
ground state of the Hamiltonian. It leads to a non-zero value of vacuum Bogoliubov 
quasiaverage of cr x which is spontaneous magnetization: 

{o**)lCY = XY<l^' y |>XY , (^)xY = V2 ( ^"y ) V8 , <^>XY = . 

The matrix element xYd^DxY & Ye spontaneous magnetization because it appears 
as 0-particle contribution from R-sector in the long-distance expansion of two-point 
correlation function (cXq<7^) in the ferromagnetic phase. 

6. Asymptotics of the two-point correlation function of the infinite 
XY-chain 

We are interested in the asymptotics of the two-point correlation function when the 
distance d between the correlating spins and is large (and, of course, the length n 
is much larger). We show how to re-derive the formula from [3] without the use of Szego 
theorem and Wiener-Hopf method for finding the asymptotics of Toeplitz determinant. 

For definiteness we consider the paramagnetic phase. In this case the intermediate 
states between two spin operators are odd-number particle states from R-sector. Let us 
show that the main contribution is due to 1-particle states (the contribution of many 
particle states is exponentially suppressed for large d). We have^] 

(*0*d) = Ns(Kb>R • r(pK\)nS 
PGR 

+ Yl NS(K|P1,P2,P3)R ■ r(Pi,P2,P3\0-}\)nS + ' ' ■ 

{Pl<P2<P3}GR 

= E eipd l N s<K»R| 2 + E e i{pl+P2+P3) %s(K\Pi,P2,P3)R\ 2 + --- , 

PGR {pi<P2<P3>GR 

where |)ns is the vacuum state. The idea is to make a transformation which corresponds 
to change of the direction of evolution to the transverse direction in the lattice 
formulation of the Ising model. Let us estimate the L-particle contribution (in the 
limit n-^oowe use integrals instead of sums): 

cosh 2 K; Yl ^ P1+ '" +PL)d \m(K\Pi, ■ ■ -,Pl)r\ 2 = t 

{pi<---<Pi}GR 

§ In this section all the eigenstates of XY-chain Hamiltonian arc labeled only by NS or R depending 
on the sector. 



t y t y 
f - f-l 



1 \ / 
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X 



_j_ r ... rfr .^., ::t - s.)W n (37) 



Now we change the variables in the integrals from the momenta {pi} to parameters {it;} 
on the elliptic curve by means of relations from Sect. 13.41 and from [32] 

. ty 1 ~ ty dpi sin^ t^ - 1 

UUl = — : ; r , , . , , 7 = ; SnlM; — Uv , K ) . 

4 sinh 7 (pO sinh lW±lW 2 k yi ' 

Thus the L-particle contribution is 

/i.-cV 2 ' 2 1 4 1 - 1 /i-'-fA 1 ' 1-1 ' /- 2 * /- 2K ., 



x n -Am + i^72) d [n \{ Ul + iir'/2) i /2 + j] \{ui + i^72)- i / 2> ) n 

Z=l \ I I J KV 



Shifting the contours of the integrations by —iK'/2 we get 

( t v -t? \ L2/2 1 4^ L ™ [™ [** 



x n (n a ( m ') i/2 + n a m~ i/2 ) n sn2 ^ - ri ) • 



Z=l \ Z Z / Kl' 

With the use of the relations X(ui) = e lpl , z{u{) = e~^ Pl \ 

t- 1 - 1 dp, sin^ V 1 - t y _ 

= : ; — r , _,_ , , _,_ , = — ; SnlW; — Uv . k ) , 

4 sinh7(p ; ) ginh 7fa)+7^) 2 k l/ ; ' 

the L-particle contribution to the two-point correlation function can be rewritten in the 
terms of energies and momenta corresponding to the transverse direction of the Ising 
lattice: 

^ ,'tx-t- x \ L2 ' 2 1 r r d PL 



ty-ty 1 / (27r) L L! J_ n J_ n sinh 7(p x ) ' ' ' sinh 7(p L ) 



x e -rfE i=1 7ta) cos 2 



U=lPlTT Sm 2 (3g) 
, sinh 2 7 ^ p ^ +7 ( p '^ 



z<z' 0±±±±± 2 



In the limit d — > oo the main contribution of these integrals is given by the neighborhood 
of zero momenta, where j(p) has a minimum. When the parameters of the model 
are non-critical we have < exp(— 7(0)) = a 2 < 1 and the asymptotics of L-particle 
contribution to the two-point correlation function is proportional to e -QLL7( - cl ). Thus the 
l-particle contribution dominates in the paramagnetic phase. 
From ( 1371) for the l-particle contribution we have 

cosh 2 K; ^ J- (^|^) 1/2 £ dp ^ coth M . 
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Now we take into account (I14p . (1151) in order to rewrite the formula in terms of a and 
b. With z = e ip we have 

7(p) a af (b 2 — z)(b 2 — z" 1 \ 1 ^ 2 



b p b \{a 2 - z){a 2 - z~ l 
We rewrite the integral contour integral: 




^J-n 2 27rijj^ =1 y (a 2 - z)(a~ 2 - z) 

To fix a branch of function under the integral we make two cuts: from b 2 to a 2 
(0 < b 2 < a 2 < 1) and from a~ 2 to 6~ 2 (1 < a~ 2 < b~ 2 ). We move the contour to 
be along the cut from b 2 to a 2 : 



7T 



6= 



{a 2 — s)(a~ 2 — s) 

Since we consider the limit d — > oo, the the main contribution to the integral is given 
by the neighborhood of a 2 due to the factor s ' -1 . We change the variable of integration 
to t = a 2 — s: 



Wo V a 2 - a 2 + t 

and use the expansion at d — > oo and £ ~ 

( a 2 _ = e (d-l)(loga 2 +log(l-t/a 2 ) ^ fl 2(d-l) e -(d-l)t/a 2 

to get the leading asymptotics by steepest descent method: 



1(d) ~ a 2(d ~ 1} (d 



. 1/2 / a 2 (a 2 -6 2 )(5- 2 -a 



2> 



ir(a 2 — a 2 ) 



Subsequent terms of the asymptotics for the 1-particle contributions are given by (4.22) 
of [3j. Finally the leading term of asymptotics of the two-point correlation function is 

M ~2{l—^) ^lq~ ad 71 ■{ ) ' (39) 

where we used (JUJ) , (J23D and fl33]). It coincides up to factor (-1)74 with (4.25) of [3] 
after identification Ai = 6~ 2 , A2 = a 2 . The origin of this discrepancy factor is the 
following: In [3] the spin operators are 1/2 of Pauli matrices. It gives factor 1/4 for the 
two-point correlation functions. Also in [3J, the coefficients at cr^a^ +1 and o~\a y k+1 in the 
Hamiltonian are positive. In this paper they are negative. To change the signs we make 
automorphism of the algebra of spin operators: <j£ — > (— l) k o~l, a\ — > (— 1) o\ —> a\. 
This map changes the signs of the mentioned coefficients and leads to the factor ( — l) d 
in the two-point correlation function. To rewrite the correlation function (1391) in terms 
of h and x one needs (TT9l). 
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Another way to calculate the two-point correlation function is to use integral 
representation for the L-particle contribution ( 138]) in terms of energies in transverse 
direction on the Ising lattice: 

cosh 2 K* ■ (a x a x d ) ~ ±- V~ ^ / dp %- , , cos 2 V - 
y ^ ° d ' 27rt-i-t y J_„ ^sinh7(p) 2 



£ t x l ~ t x r J„ „-<i(7(0)+7"(0)p 2 /2) ^ J_ t ^ ~ tx e / 271 " 



dp e 



" 2tt t^ 1 - sinh7(0) J_ n y " 2tt t~ l - t y d 1 / 2 sinh7(0) ]] f'(0) ' 

We use also ( ITS"]) and 

7"( ) = 1 jLz1^ . 1 , e -?(o) = a 2 S inh7(0) = ^ ~ - 
' V ; V ^ sinh7(0) ' [V ; 2 

to obtain by steepest descent method the same asymptotics of the two-point correlation 
function ( 1391) . 

To obtain more precise asymptotics, we change the variable of integration s = 
exp(— 7(p)), —ds/s = dj(p) and use the dispersion relation (!T8~]) to get 

dp d7 V 1 ~ ^ _ rfs V 1 ~ 



sinh7(p) ship t x 1 — t x s ship t x 1 —t a 



It gives 



,2 



£-i _ ^ r n e -i{p)d PI" p 
dp cos 2 — = ds 3 d " 1 cot - 



t y l - t yJ-^ sinh7(p) 2 J fe2 2 
From (IT51) we have 



p (s — b 2 )(b~ 

cot 



2 y (a 2 - s)(a~ 2 - s) ' 

Substituting this expression into the integral we get the formula coinciding with formula 
(4.21) from [3]. It can be expanded as described there to obtain the subsequent terms 
of the asymptotics. 



7. Discussion 



Using the results of pTJ on the separation of variables method for the Baxter-Bazhanov- 
Stroganov model, in this paper we derived matrix elements (form-factors) of spin 
operators cr x and er y between the eigenvectors of the Hamiltonian ([1]) of finite (length 
n) XY-chain. The final formulas are ( 13TI) . ( 132]) and ( l20i) . In the limit of infinite 
chain the formulas are reduced to (|34|) . (1331) and (133]) and allow to re-derive the 
asymptotics of correlation function (er x er x ) at d — > oo in a simple way. It is interesting to 
derive analytically the long-time and long-distance asymptotic behavior of the a x spin 
correlations at finite temperature observed recently in [38] by numeric analysis. 

Usually the standard efficient way to study the Ising model and XY-chain is to use 
the language of fermions. Recently the factorized formulas for the matrix elements of 
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spin operators were re-derived for the quantum Ising chain [53], for the Ising model on 
two-dimensional lattice [39], for general free fermion model on two-dimensional lattice 
and quantum XY-chain [29] using algebra of fermion operators. 
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